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Scanning tunneling spectroscopy of superconducting LiFeAs single crystals: Evidence 
for two nodeless energy gaps and coupling to a bosonic mode 
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The superconducting compound, LiFeAs, is studied by scanning tunneling microscopy and 
spectroscopy. A gap map of the unreconstructed surface indicates a high degree of homogeneity 
in this system. Spectra at 2 K show two nodeless superconducting gaps with Ai = 5.3 ±0.1 meV 
and A2 = 2.5 ± 0.2 meV. The gaps close as the temperature is increased to the bulk T c indicating 
that the surface accurately represents the bulk. A dip-hump structure is observed below T c with an 
energy scale consistent with a magnetic resonance recently reported by inelastic neutron scattering. 

PACS numbers: 74.55,+v, 74.20.Mn, 74.20.Rp, 74.70.Xa 



In the new family of high temperature superconduc- 
tors, the iron pnictides, a general consensus is emerging 
in favor of an s± symmetry [l|-[3[ although other paring 
states such asp-wave 0, [3 and d-wave have also been 
suggested [3, Q. A challenge in unambiguously identi- 
fying the pairing state in this class of materials is that 
experimental investigations are occurring against a back- 
drop of variations in sample purity and quality. In partic- 
ular many of these compounds have intrinsic limitations 
due to high defect density from cation doping in the bulk, 
or structural or electronic reconstructions which compli- 
cate surface sensitive investigations 0, 0] • 

A particularly interesting compound among the pnic- 
tides is LiFeAs which is superconducting without cation 
substitution [!, This potentially places it in the 
same position that YBa2Cus07_^ (YBCO) holds in the 
cuprates |Io|, a stoichiometric superconductor that can 
be chemically and structurally perfect enough to avoid 
artifacts arising from disorder. LiFeAs has the addi- 
tional advantage of a natural cleaving plane, exposing a 
non-polar surface that does not undergo reconstruction 
making it well suited to surface sensitive spectro- 
scopic studies such as angle resolved photoemission spec- 
troscopy (ARPES) [HI, [l3|| and scanning tunneling mi- 
croscopy and spectroscopy (STM and STS) much 
like the cuprate Bi 2 Sr 2 CaCu208 + x (BSCCO) [13 EH. 



In this letter, we show through STM that the sur- 
face is clean and unreconstructed. Through spatially re- 
solved STS measurements, we find that the gap structure 
is extremely homogeneous, presenting an opportunity to 
study a clean and well-defined system. STS acquired at 
2K reveal two nodeless gaps, consistent with a multiband 
s± pairing state [l|. We find that the temperature depen- 
dence of the gap is BCS-like, in contrast to the fluctuation 
driven transition [19l| a nd pseudogap phase present above 
T c in the cuprates [20|. All of these simplifying features 
offer a system in which to study a feature LiFeAs does 



have in common with the cuprates; a pronounced struc- 
ture above the superconducting gap, indicating strong 
coupling to boson modes. 
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FIG. 1. (a) 80 x 80 nm 2 scanning tunneling topographic image 
of LiFeAs (V B = 50 mV, I T = 20 pA) at T=4.2K and (b) 
6.8 x 6.8 nm 2 atomic resolution topographic image of LiFeAs 
(V B = 40 mV, I T = 100 pA). (c) Schematic of the crystal 
structure from the cleaved (001) direction and an edge-on 
view, (d) dl/dV spectrum in a range from -550 mV to +550 
mV. 

Single crystals of LiFeAs were grown by a self- flux tech- 
nique. LiaAs and FeAs, pre-synthesized from Li (99.9%), 
Fe (99.995%) and As (99.9999%), were mixed in a com- 
position of 1:2 and sealed under 0.3 atm Ar. The mix- 
ture was heated to 1323 K for 10 hours, then cooled to 
1073 K at 4.5 K/hour. Finally, the samples were addi- 
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tionally annealed at 673 K for 12 hours before being re- 
moved from the furnace. Single crystals with dimensions 
up to 2x2x0.2 mm 3 were mechanically extracted from 
the flux. Lattice parameters a = (3.777 ± 0.004)A and 
c = (6.358 ± 0.001) A were determined by x-ray diffrac- 
tion, and T° nset = 17 K with a transition width of IK 
was determined by SQUID magnetometrv (see Fig. [3(c)), 
consistent with stoichiometric LiFeAs [9|. 

All measurements were performed in an ultrahigh vac- 
uum liquid helium temperature STM. The electrochemi- 
cally etched tungsten tips were Ar + sputtered and ther- 
mally annealed at the beginning of each set of experi- 
ments. The samples were cleaved in-situ at a tempera- 
ture of ~20 K. Reproducible datasets were obtained from 
studies of four different samples, some of which have been 
recleaved and studied with individually prepared tips. 
All spectra shown here were acquired by numerical dif- 
ferentiation of the I-V sweep. 

Cleavage most likely occurs between two weakly 
bonded lithium layers [11], resulting in a surface that con- 
sists of a lithium layer, with arsenic then iron beneath, 
shown in Fig. H(c). A topographic image of an 80 x 80 
nm 2 area is shown in Fig. H(a), in which several types 
of defects are visible. The defect density measured from 
STM images is 0.0020 ± 0.0005 per LiFeAs formula unit. 
Such low defect density further confirms the stoichiomet- 
ric composition of our LiFeAs crystals. Figure QJb) shows 
an atomic resolution topographic image with a size of 
6.8 x 6.8 nm 2 (V B = 40 mV, I T = 100 pA). The inter- 
atomic spacing of the square lattice is (3.74 ± 0.03) A, 
which is consistent with the periodicity of either lithium 
or arsenic at the surface (see Fig. 0Jc)), but not with the 
iron lattice. 

Figure H(d) shows a dl/dV spectrum in a bias range 
from —550 mV to 550 mV with a superconducting gap 
at Ep. The Fermi energy is located in a region of low 
local density of states (LDOS) enclosed by enhancements 
below about —250 meV and above 10 meV. The sharp in- 
crease of the LDOS just above the Fermi energy reaches a 
plateau at 35 meV. Similar features have been explained 
by a surface state, as summarized in Q. However, a 
surface state in LiFeAs is unlikely according to density 
functional theory 
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H2||. The overall U or V shaped 
background DOS is a generic feature of iron-based super- 
conductors, that has been attributed to their semimetal- 
lic character caused by a small overlap of the bottom of 
electron bands and the top of hole bands close to Ep 

To gain deeper insight into the spatial homogeneity of 
the superconducting gap magnitude and in order to study 
how defects influence the gap magnitude, we recorded 
400 x 400 STS spectra within a bias range of ±15 mV over 
an area of 90 x 90 nm 2 (see the topography in Fig. [5(a)). 
Figure [2(b) shows the 150 measured STS spectra along 
the 33 nm white line indicated in Fig. [5(a), demonstrat- 
ing a striking homogeneity of the superconducting gap 
in defect free areas. A gap map from the whole image is 
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FIG. 2. (a) 90 x 90 nm 2 topographic image (Vb — 15 mV, 
It = 100 pA) of LiFeAs. (b) 150 individual dl/dV spectra 
(T=4.2K) from the defect-free 33 nm long line marked in (a), 
(c) Gap map of the same region in (a). A hpp corresponds to 
half the energy separation between coherence peaks. Both, 
the gap map and the histogram of A hpp shown in (d) reveal a 
high degree in homogeneity of the superconducting properties 
in LiFeAs. 



extracted from A hpp , which is half the energy separation 
between apparent coherence peak maxima. A hpp gener- 
ally overestimates the magnitude of nodeless gaps due to 
thermal broadening (here T=4.2 K). However, A hpp does 
allow for an analysis of the spatial variation, independent 
of any specific model. The gap map Fig. [5(c) and the 
corresponding histogram Fig. [5(d) indicate a consistent 

magnitude with a mean value A =6.07 meV and a 
small standard deviation a = 0.08 meV. Much of the 
variation comes from spectra near defects, where A hpp 
is reduced to a minimum of 5.7 meV. The homogene- 
ity of the superconducting gap measured in LiFeAs by 
STS with <ta/A ~ 1.3% is in contrast to most other 
high temperature superconductors. For comparison, the 
122 iron arsenide compound BaFe1.8Coo.2As2 was found 
to have a less homogeneous gap with oa/A ~ 12% 
2l|. An extreme example is the much-studied cuprate, 
BSCCO [22(| .where local inhomogeneities in doping re- 
sult in strong variation in the gap magnitude. Hence, 
free from surface reconstruction or disorder induced by 
local dopants, LiFeAs presents an ideal system for surface 
sensitive spectroscopic investigations. 

Taking a closer look at the STS in the vicinity of the 
superconducting gap in Fig. [3(a), two nodeless gaps are 
clearly resolved at 2 K. The two-gap superconducting 
excitation spectrum can be fitted within the framework 
of BCS, where the normalized quasiparticle density of 
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FIG. 3. (a) Two isotropic gaps fit (blue line) on top of mea- 
sured dl/dV spectra (black line, 2 K). The dashed lines indi- 
cate the fitting bias range from -6.8 mV to 6.8 mV. (b) Tem- 
perature dependence of the superconducting gap measured 
by STS between 2 and 17 K. (c) The large gap determined 
by isotropic s-wave fits (red error bars) generally follows the 
temperature dependence predicted by the BCS weak coupling 
limit (dashed black line) . The development of the smaller gap 
(blue error bars) is obscured by thermal broadening at ele- 
vated temperatures. The bulk Tc of 17 K probed by SQUID 
magnetometry with a 1 G magnetic field (gray dots and right 
y-axis) agrees with the surface critical temperature, demon- 
strating agreement between surface and bulk properties. 



states N(E) of a superconductor is defined as [2S 

N s (E(k)) 



N(E) = J2- 



E(k) 2 - A(k) 2 



(1) 



with the superconducting DOS N s and normal states 
DOS 7V n being a function of the Cooper pair energy E(k) 



and the single particle energy y E(k) 2 — A(/c) 2 , respec- 
tively. If measured by the tunneling method, this can be 
described by Dynes' formula [24j . 



N(eV) = Ylf Re 



x (E - iT) x 



df(E-eV) 
dE 



i J(E-iT)i-\A t (kW 



dE 



(2) 



where f(E — eV) is the Fermi-Dirac distribution function 
and Wi is the weight of the contribution from the i th 
with the constraint of w\ 



gap 

w 2 = 1 25]. The effective 
damping term T(E) = aE is used to properly represent 
zero DOS at E F [26]. We fit the 2 K STS over a bias 
range from -6.8 mV to 6.8 mV under the assumption 
that the normal density of states is linear over the small 
energy range examined: N n {E) = a + b x E. 

Two different A(fc) superconducting gap models were 



considered. One consists of two isotropic gaps, yield- 
ing Ax = 5.33 ± 0.10 meV, A 2 = 2.50 ± 0.15 meV, 
k;i=0.89 and energy dependent F = 0.13 x E. The 
other consists of two anisotropic gap s with four fold 
symmetry as observed by ARPES [13|, yielding Ai = 
5.33 x (1 + 0.09 x cos(40))±O.l meV (Af ax = 5.8 meV), 
A 2 = 2.50 x (1 + 0.20 x cos(4(9 / )) ± 0.20 meV, wi=0.87 
and r = 0.10 x E. Both gap models fit the 2 K STS 
very well within the fitting range (see for example the 
two-isotropic-gap model fit in Fig. l2fb)) and give gaps 
that are consistent with recently reported values obtained 
from STS 15, 16|, giving reduced gaps of =7.3 and 
fc 2 ^ = 3.4. However, both of the fits clearly fail to rep- 
resent the measured spectra outside the fitting range due 
to additional above gap features, which will be discussed 
further below. 

The gap magnitudes obtained from fitting agree well 
with the largest and the smallest of the four gaps de- 
termined by ARPES [13]. Tunneling into the two elec- 
tron pockets, containing the other two gaps, is expected 
to be strongly suppressed because of the larger in-plane 



momentum |fc||| [6|, |27j. The small anisotropy factors ob- 
tained from the fit to two anisotropic gaps also agree well 
with ARPES results [Hj, reinforcing the consistency of 
the two surface sensitive measurements in LiFeAs. How- 
ever, unlike the gap shapes previously reported by STS 
[lij . the gaps shown in Fig. [3ja) and (b) are fully open 
and symmetric even at elevated temperatures. 

Fig. [3]^b) shows the temperature dependence of the 
STS spectra from 2 K to 20 K in the same defect-free 
region. Each spectrum is the average of 16 spectra ac- 
quired from a 4 x 4 nm 2 area. The superconducting gaps 
are visible up to 16.5 K and disappear at 17 K, mani- 
festing the same T c seen in susceptibility. Based on the 
simpler isotropic gap assumption, which adequately rep- 
resents the gap size given the small anisotropy factors 
found, the temperature dependence of the gap ampli- 
tudes was extracted, and is plotted with the Meissner 
transition measured by SQUID in Fig. [3fc). The tem- 
perature dependence of Ai follows BCS theory and the 
gap closes at the bulk T c value. 

These results and the apparent absence of a surface 
state or electronic reconstruction 11, [l2| suggest that 
the surface behavior echoes the bulk. 

To more carefully examine the structure surrounding 
the superconducting gap, the STS below T c are normal- 
ized by the normal state conductance, dl /dV(V)\uK 
according to Eq. [TJ shown in Fig. H^a), and the sec- 
ond derivative was calculated numerically, shown in Fig 
2t[b). Here, features higher in energy than the coherence 
peaks are clearly visible, diminishing as the tempera- 
ture approaches T c . These features consist of a relatively 
well-defined dip below the normal state conductance, fol- 
lowed by a broad hump. These can be characterized by 
three different energies: Ejj ^10meV, Ej ^12meV, and 
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or shoulder. Regardless, features caused by coupling to 
a bosonic mode are expected to appear at the mode en- 
ergy shifted by the gap (CI 



42|. Given the dif- 
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FIG. 4. a) Normalized STS in a temperature range between 
2 and 16.5 K. A dip-hump feature is clearly visible at 2 K 
but decays with increasing temperature, b) The normalized 
second derivative of I(V) at 2 K reveals the position of the 
main features. The dip in the normalized dl/dV is at Ed ~ 
±10 mV and the inflection point/kink in dl/dV is at Ei « 
±12 mV (dashed line). 



Eh ~15-20meV, the energy of the dip, inflection point 
between dip and hump, and the hump respectively. 

While the apparent bias asymmetry in the coher- 
ence peaks and dip-hump features are reminiscent of the 
this may also arise from the normalization 



cuprates 

to the 17K spectrum which exhibits a hump at around - 
5meV and a steep rise at positive bias, both of which are 
thermally broadened compared to the low temperature 
spectra. 

The features observed in LiFeAs, characterized by a 
dip at rsj 2 A from Ep followed by a broader hump, bear 
stri king resemblance to those observed in the cuprates 
17l l28-[30|. In past studies on superconducting cuprate 
materials, these features have been attributed to several 
different origins due to the large variety of competing 
effects at similar energy scales. These include inelas- 



tic tunneling effects [3l|, |32| , band structure effects [33| , 
though most attribute these spectra features to either 
pairing [28|, H, l34] - [37| or non-pairing boson interactions 
[38| . Other explanations are based on the pseudogap ob- 
served in the cuprates j39l Eoj , but these are likely absent 
in the iron arsenides [6| . In our data, the reduction of 
the tunneling conductance below the normal state indi- 
cates that inelastic tunneling effects are not responsible 
for these features, and band structure effects are unlikely 
in our spectra since the features disappear above T c . 

Thus, we turn our attention to possible boson inter- 
actions. In the well-established framework of phonon- 
based pairing in an s-wave superconductor [4l|, the en- 
ergy dependence of the gap leads to an initial peak or 
shoulder outside the quasiparticle coherence peaks due 
to increased pairing strength, followed by a dip as the 
interaction switches from attractive to repulsive at the 
mode energy. Our spectra do not show this initial peak 



A) |2a 

ferences between our spectra and the classic phonon cou- 
pling case, combined with the lack of strong features in 
the phonon spectrum of LiFeAs below ~14meV 43, 44j . 
it seems unlikely that the features we observe arise from 
phonon mediated pairing. Spin fluctuation mediated su- 
perconductivity has been suggested in the pnictides [l[, 
a theory also supported by STS data of SmFeAsOi-^F^ 
(45| . Indeed, recent reports from neutron scattering have 
indicated a broad magnetic excitation peaked around 5- 
8meV (isl, |4?J , corresponding well with the energy of the 
dip shifted by the large gap, Ejj — Ai ^5meV. 

The cuprates and the arsenides thus share a couple 
of features that suggest a common origin for the dip- 
hump. Both have a mode developing in the spin fluctu- 
ation spectrum when the superconducting gap appears 
below T c 0, HqI • Both also have considerable damping 
of this mode due to non-zero density of states; in the 
cuprates this is due to nodes in the d-wave gap and in 
the arsenides this can arise from anisotropy as well as 
a finite density of states once above the energy of the 
small gaps associated with the multiband nature of the 
s± state. This may reasonably result in similar spectral 
shapes, however unambiguous assignment of the spec- 
tral features will require a proper microscopic theory in- 
cluding the pairing symmetry of LiFeAs. The energy 
scale of the features found here also draw a parallel to 
the cuprates. Zasadzinski et al. showed a clear propor- 
tionality of the boson mode energy Ct with T c , that also 
agreed in magnitude with the magnetic resonance from 
inelastic neutron scattering for a wide range of T c from 
overdoped to underdoped BSCCO[29|. They found that 
Ct/ A « 1 for optimally doped samples and generally that 
Ct w bksTc, confirmed for a wide range of cuprate super- 
conductors j48|. The feature observed here for LiFeAs 
similarly sits at Ct/A « 1 and lies close to the 5fc£T c = 7 
meV, indicating a similar underlying mechanism. 

The characteristics of LiFeAs presented here demon- 
strate that this material provides a comparatively sim- 
ple system in which to study high-T c superconductiv- 
ity. In stark contrast to the best cuprate materials, the 
superconducting gap remains remarkably homogeneous 
over large areas in LiFeAs. Additionally, the presenta- 
tion of a non-polar cleaved surface, that does not re- 
construct and accurately represents the bulk properties, 
makes it ideal for surface sensitive studies. Although this 
material shows multiple superconducting gaps, they are 
without nodes, and exhibit BCS-like temperature depen- 
dence. Yet, this material shows all the signs of strong 
coupling with a relatively large reduced gap of 7.3, and 
strong above gap features corresponding closely in energy 
with a magnetic resonance recently reported. The sum 
of these features offers a clearer view into the quasiparti- 
cle phenomenology and should serve as a testbed for new 
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understanding of superconductivity. 
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